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I.      INTRODUCTION 

The  possibility  of  making  aocurat©  mechanical  forecasts   of  weather 
phenoinena  by  straightforward  solutions  of  the  basio  equations   governing 
atmospheric  motions  has   long  intrigued  meteorologists.     Theoretically, 
such  a   solution  to  the  problem  is  possible  but  there  are  so  many  varied 
and  complex  forces  at  work   in  the  atmosphere  that  all  attempts  so  far  to 
effect  a   complete   solution  to  these  equations  have  met  with  failure* 

One  of  the  more  important  difficulties  in  the  way  of  reaching  satis- 
factory solutions  has  been  the  non-linearity  of  the  equations  and  no 
standard  methods  are  available  to  integrate  nonlinear  partial  differential 
equations •     Therefore  the   introduction   of  the   so-called  ''perturbation" 
method  was  made  which  consists  essentially  in   linearizing  the  equations  by 
assuming  meteorological  quantities  having  basis  values  unchanging  with  time, 
plus   small  perturbation  values  whose  second  order  terms  can  be  neglected* 

One  approach  to  quantitative  forecasting  has  been  the   linearization 
and   solution  of  the  absolute  vortloity  equation.     This  was   introduced  by 
Rossby     flO]    ,  who  obtained  a  solution  for   the  motion  of  sinusoidal  waves 
of  infinite   lateral  extent  in  a  horizontal  plane  and  found  that  the  velocity 
of  propagation  of  such  waves   is   given  by  the  well  known  formula 

where  C  is  the  velocity  of  the  waves  toward  the  East,!/  is  average  speed  of 
the  westerly  current,  L.  is  wave   length  and  ^  is  the  meridional  rate   of  change 
of  the   Coriolis  parameter*     A  considerable   literature  has  grown  out  of  the 
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interpretation  of  this  formula  on  the  classification  of  large  scale 
atmospherio  systems   in  terms  of  the  so-called  tonal  index«  which  is  a 
measure  of  the  strength  of  the  zonal  current* 

Hauri»its  extended  Ro8sby*s  method  and  obtained  solutions  for  the 
motion  of  waves  of  finite   lateral  extent  on  a  horitontal  plane    £  5J    and 
on  a   sphere     [6j   .     Craig     [  l]     obtained  solutions   of  the  vortioity  equa- 
tions  in  complete  form  without  linearisation  for  the  plane  and  the  sphere. 
These  solutions  differed  from  those  of  Eaurwitz  by  the  absence  of  a  zonal 
velocity  term.     Neamtan     [  7j     then  made  a  treatment  of  the  vortioity  in 
the  same  manner  as  Craig,  viz  by  the  use  of  stream  functions  and  obtained 
identical  solutions  to  those  of  Haurwitz*     He  showed  that  Craig* s  solutions 
were  incomplete ,   causing  absence  of  the  zonal  index  term* 

Forsythe    [  3}     developed  formulas  for  the  speed  of  propagation  of  waves 
with  changing  shape  by  use  of  the  scalar  relative  vortioity  as  an  identifiable 
property.     However,   he  did  not  make  any  test  on  his  formulas  to  determine 
their  usefulness   in  practical  forecasting* 

Rossby  and  co-workers     [11]     supplemented  his   initial  formula  based  on 
the  conservation  of  absolute  vortioity  with  a  more  general  technique  which 
would  be  applicable  to  arbitrary  initial  streamline  pattern  so  that  the  dif- 
ficulty of  defining  a  prevailing  wave  length  would  be  circumvented* 

All  of  the  above  solutions   of  the  vortioity  equations  were  predicated  on 
the  assumption  of  friotionless,  n on -diver gent,   and  autobarotropic  flow. 
Further,    in  order  to  obtain  a  solution  for  relative  vortioity  which  could  be 
useful  in  forecasting  trajectories   it  was  necessary  for  Rossby    Cll]     to  as- 
sume no  horizontal  shear  of  the  wind  and  stationary  pressure  systems* 
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In  thie  treataient,  a  general  solution  of  the  non-linsar  differential 
equation  of  the  vertical  component  of  absolute  vorticity  iwill  be  obtained 
without  assuming  non-divergent  flow  or  baro tropic  conditions.  A  solution 
is  obtained  for  (1)  a  surface  of  constant  height,  and  (2)  a  constant  pres- 
sure surface*  A  discussion  of  mean-value  constants  for  the  wind  shear  term 
and  for  the  effect  of  moving  pressure  systems,  by  means  of  which  the  vorticity 
equation  could  be  integrated  without  disregarding  these  two  terms,  is  included* 
A  chapter  is  devoted  to  examining  this  solution  with  regard  to  its  adaptation 
to  forecasting  of  upper-air  flow  pattern  by  the  trajectory  method.  In  addition, 
special  forms  of  the  vorticity  equation  resulting  from  various  assumptions  are 
developed  and  discussed* 

Gradient  flow  is  assumed  throughout  this  investigation.  The  use  of 
gradient  flow  permits  a  special  handling  of  the  divergence  term.  The  wind 
of  course  is  not  always  gradient  even  in  the  free  atmosphere  but  is  closely 
approximated  by  the  gradient  wind  at  elevations  greater  than  1000  meters 
above  the  ground.  However,  under  conditions  of  rapidly  changing  pressure 
gradient  this  close  relationship  between  observed  wind  and  gradient  wind 
is  greatly  modified,  and  is  due  to  the  fact  that  the  motion  is  not  under 
balanced  forces*  Under  these  conditions  there  is  a  velocity  component 
along  the  isallobaric  gradient.  Consequently  the  use  of  gradient  flow  in 
the  vorticity  equation  would  generally  be  in  error* 
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II.   SOLUTION  OF  THE  VORTICITT  EQUATION  FOR  GRADIENT  FLCW 

In  obtaining  solutions  to  the  differential  equation  of  the  vertical 
component  of  absolute  vortioity  as  few  assumptions  as  possible  will  be 
employed.  Initially,  only  the  following  assumptions  are  madet 

a.  Friction  is  negligible 

b.  Vertical  relooitles  are  negligible. 

Other  restrictions  will  be  added  later  to  obtain  solutions  of  particular 
oases.  The  general  equation  of  vortioity  to  be  solved,  which  is  due  to 
Bjerknes,  is: 

^= -:r«'7.\/4  7ivr.  \<5,     .  (2.1) 

It  is  noted  that  this  equation  consists  of  two  terms,  a  term  which 
represents  the  horizontal  divergence  of  the  gradi«rit  wind  and  a  sole- 
noidal  term  which  expresses  the  component  of  the  geostrophic  wind  along 
the  gradient  of  temperature* 

Since      N/^,  =  -  -^  ^-px  ^  ^ 

^._  :^^v_vr.o<f^-^xi^j        •  .         (2.2) 

According  to  Taylor,  TlSj   the  following  expression  for   V./?\/ can 

A  ^' 

be  derived: 

which,  solved  for  V-V  ,  yields 

-         :^  ^    • 
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Now  if  we  put  this  value  of  V  V  in  equation  (2*2)  and  divide  by  ^ 
we  obtains  :^  <>-?--  /'-»^'"^/' S  V^  7.  ' 


^.    A--*'^    /'■'      '^    ^  '       •  '       '>    / 


Now  v^-  vf^  VT  , 

80  that 

:f.  ;F--  +  ~>      r    +-^  T     -^"-   .(2.4) 

\/.^Ps.o     oeoause         raotors/,  perpendicular  •     Henoe 


-/ 


^  ^^        -     -r  ^  -r         -^ 

Now  if  we  make  the  substitution 


we  obtain 


Rearranging  and  collecting  terms 
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Integrating  along  the  path  of  the  particle. 
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and  since      V^  =   ^-  > 


(2.6) 


NofwJ^»/)f;(/y4,£:,  where    ^     is  the  horizontal  shear  of  the  gradient 
wind,  and  the  relation  of  ourvat\we  of  streamlines  to  curvature  of  trajec- 
tories  i^  J^\/^  kiv-f^Jt-    •     Hence  equation  (2»6)  can  be  written,  after 
integrating  and  rearranging  terms. 


where  the  subscript  0  indicates  the  initial  values  of  the  variable   quantities* 

We  note  from  equation   (2.5)  that  on  following  a  particle  along  a  path 
from  an  initial  point  to  some  other  point  on  the  path,   the  ratio  of  their 
absolute  vortioities  is  equal  to  the  ratio  of  their  temperatures  times  the 
inverse  ratio  of  terms   involving  curvature  and  shear  times  an  exponential 
term. 

Equation  (2«6)  was  derived  without  assuming  either   zero  divergence  or 
a  barotropic  atmosphere.     Furthermore,  no  restrictions  have  yet  been  placed 
on  the  shear  of  the  wind  or  changing  streamlines.  - — 

Ho,  if  we  Ut     4,^m-(ki-l/,)v]         i^     I,     %,7if^'l^-s)\/ 
.and      ^  -    mean  value  of   ^     along  the  path,   then  '^/o 
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If  we  assume  a  stationary  pressure  system,  t^x"  -  ^   »  **'  ^®  most 
investigations  of  vortioity,  then  equation  (2 06)  becomes 


2>/^  •  (2.8) 


1.     Discussion  of  Orders  of  Magnitude  of  the  Mean  Value  Constant  in 
Equation  (2*8)  for  Stationary  Systems. 

If  we  compare  the  relative  order  of  magnitude  of  the  terms  contained  in 
the  exponential  part  of  equation   (2*8)   some  conclusions  can  be  reached  re- 
garding the  value  of  the  term  -^ ^      so  that  a  constant  can  be  assigned  to 
this  term  which  will  represent  its  average  value  over  most  mid-latitude 
paths •    With  such  a  constant  the  exponential  terms  can  then  be   integrated* 

An  examination  of  700  mb.  charts  appears  to  demonstrate  a  radius  of 
curvature  of  most  waves  of  the  order  of  600  miles*     If  we  take  an  average 
value  of  20  knots  for  the  wind  velocity  the  order  of  magnitude  for  the 
curvature  term  will  be  10""^  sec"^«     Estimating  an  order  of  magnitude  for 
the  shear  term  in  a  similar  manner  we  use  for  an  average  shear  a  change  of 
20  knots  in  a  distance  of  600  miles*     This  gives  an  order  of  magnitude  of 
10"°  seo"^  for  the  shear  term*     An  average  value  for  the  coriolis  force 
term   A   would  be  the  value  at  45^  latitude*     At  this  latitude  the  coriolis 
term  is  almost  exactly  equal  to  10"^  sec"^.   .Therefore,   the  order   of  magni- 
tude for  the  curvature  and  shear  terms  will  closely  approximate   10"°  sec" 
and  for  the  coriolis  term  will  be  10"^  seo"l*     It  is  seen  from  this  qiali- 
tative  comparison  that  the  coriolis  force  is  about  ten  times  greater  than 
the  other  two  terms   in  influencing  changes  in  the  absolute  vortioity* 
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llVh«n  these  value*  are  placed  in  the  expression  ZT   ^  we  obtain 

the  value  of  -'/a    when  the  shear  term  is  negative  and  the  value  of  /// 

when  the  shear  is  positive* 

Therefore  the  general  solution  (2*8)  can  now  be  written  as 

where  \^=  vf        or  - '//   • 

The  f^j       term  will  then  be  /S)        or  (J-J.  (Z/Vj  differs 

from  2»0  by  about  3%*     Now,  since  the  ratio  of  /CZT)  will  usually  be  of 
the  order  of  imity  the  error  involved  in  using  simply  /Xj    instead  of 
OtoJ      will  only  be  about  1%* 

Hence,  by  neglecting  the  exponential  term  represented  by  V   «  equation 
(2«l>)  can  be  written  simply  as 

2*   Discussion  of  Constants  in  the  Equation  for  Moving  Pressure  Systems* 
All  treatments  of  the  vorticity  equation  appearing  in  the  literature 
so  far  have  ignored  the  effects  of  difference  between  the  eurvatures  of 
streamlines  euid  trajectories  £li^  ,  and  have  assumed  a  steady  state  con- 
dition* This  implies  stationary  pressure  systems  which  of  course  is  not 
the  prevailing  case  in  nature.  The  difficulty  of  handling  ^   mathe- 
aatically  has  been  responsible  for  neglecting  this  term* 
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Perhaps  some  oonolusions  oan  be  drawn,  howoTer,  from  an  examination 
of  the  range  of  variation  of  V(^-X^/   so  that  a  mean  value  constant 
can  be  assigned  to  €     with  the  result  that  equation  (2*7)  can  be  inte- 
grated* 

From  the  relationship  of  trajectory, curvature  and  streamline  curvature 
according  to  Petterssen  ^8,  p.  225j  ^  /r/=-^(/-  ^  C^^'*')  , 
we  can  solve  for  y(kt-k^J      so  that 

Consider  now  a  wave  shaped  system  which  is  more  or  less  symmetrical 
about  a  latitude  circle  and  moving  eastward  with  speed  C.  Along  the  trough 
which  is  to  the  right  of  the  path  of  the  system  Cos  ^^^     and  is  unity  at 
bottom  of  trough*  Along  the  ridge  where  /cs    is  antioyclonio  to  left  of 
path  Cos  </*<©.   If  we  choose  average  values  of  "to.?  for  Cos  </^  ,  the  sign 
depending  on  whether  along  a  trough  or  ridge,  and  take  average  values  of 
1000  km  and  20  kts  for  /^  and  C  respectively  then  VCJt*-^)   ^"^  /^   order 
of  magnitude* 

From  previous  considerations  of  orders  of  magnitude  for  the  shear  and 
absolute  vortioity  we  arrived  at  values  of  10"^  seo"^  and  10"^  seo'^ 
respectively  for  two  terms.  Therefore  the  final  order  of  magnitude  for 
will  be  two  possible  values: 

If  the  shear  is  negative,  then  ^ r^—o2      ,  and  if  shear  is  positive. 
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Fioally,  substituting  these  two  possible  values  in  equation  (2*7) 
we  obtain 

4.    \nlV^' 


and        --.,,.   i-        ^    ^ 


Therefore,  under  the  conditions  stated  above  the  existence  of  positive 
or  cyclonic  wind  shear  cancels  out  the  effect  of  the  ^r         term  and  with 
anticyclonic  or  negative  shear  the  correction  term  6   is  significant*  The 
opposite  effects  would  occur  in  the  case  of  retrograde  waves  or  in  those 
parts  of  closed  pressure  systems  where  the  wind  blows  in  a  direction 
opposite  to  the  direction  of  motion  of  the  system* 
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III.      SPECIAL  CASES  OF  THE  GENERAL  SOLUTION  OF  THE 
VORTICITY  EQUATION  FOS  STATIONARY  SYSTQfiS 


Certain  simplifying  assumptions   can  new  be  made  with  regard  to  the 
various  terms   of  the  general  solution   (2«8)«     The  special  oases  whioh  will 
be  discussed  here  aret 

1*     Geostrophio  Flow 

2*     Wind  Flow  Parallel  to  Isotherms 

3*     No  Horizontal  Wind  Shear 

4.     Combinations  of  the  Above  Three  Cases* 
1*     Case  of  Geostrophio  Flow* 

Geostrophio  flow  prevails  where  there  exist     straight  isobars  or  contours, 
hence  the  curvature  term  becomes  zero  and  equation  (2«8)  becomes  $ 

•S,    '    X.'    ■?  ■  (3-1) 

However,  as  demonstrated  in  Chapter  II,  the  exponential  term  can  be 
neglected  with  only  a  small  error  involved  so  that  for  the  geostrophio  flow 
case  equation  (3*1)  can  be  written  as 

^.     -r^-ji        '  (3.2) 

Equation  (3*2)  can  be  written  as  a  quadratic  and  solved  by  the  quadratic 
formula  so  that 

(3.3) 
where  Cr(7)4|^\7).   . 
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Equation  (3.3)  can  now  be  expanded  binomially: 
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All  terms  beyond  the  first  in  this  binomial  expansion  are  of  the  order 
of  10"®  or  smaller  and  which,  because  of  a  muoh  smaller  order  of  magnitude 
than  for  the  shear  term,  oan  be  negleoted  so  that  for  the  case  of  geostrophic 
flow  the  absolute  vorticity  can  be  expressed  simply  by  -J^  -  V:i  ^ -^  %.^ >  ^     • 

Inasmuch  as  the  absolute  vorticity  cannot  be  sero  except  for  the  very 
improbable  case  of  geostrophic  flow  with  no  shear  and  at  the  equator,  the 
choice  of  sign  on  the  second  term  in  (3*3)  must  be  positive* 

For  practical  application  on  the  weather  chart  in  determing  how  such 
flow  would  be  affected  by  various  types  of  movement,  the  relative  vorticity 
is  examined.  For  geostrophic  flow,  (3.2)  can  be  written 


j'^frj2!:(^^^l-^  . 


A  particle  moving  northward  would  therefore  undergo  a  decrease  in 
relative  vorticity  or  be  turned  anticyolonically  even  though  there  may  be 
a  slight  positive  contribution  from  the  temperature  term.  Southward  flow 
would  produce  the  opposite  effect. 

A  particle  moving  along  a  parallel  of  latitude  would  be  influenced  only 
by  the  change  in  temperature.   If  moving  toward  higher  temperature,  there 
will  be  an  increase  in  relative  vorticity  and  particle  will  curve  to  the 
left. 
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2*  Case  of  Wind  Flow  Parallel  to  Isotherms. 

The  wind  is  blowing  parallel  to  the  isotherms  In  this  case,  which 
implies  a  barotropic  atmosphere*  Suoh  oonditions  prevail  when  systems 
are  thermally  symmetrioal*  In  this  oase  equation  (2.8)  becomes 

^*    (>K/-f;)9  (3.4) 

and  the  exponential  term  becomes  unity*  Now,  on  following  a  particle  from 
initial  point  or  path  to  some  other  point  the  absolute  vorticity  changes 
as  the  inverse  ratio  of  their  respective  (/s/-//)y  terms • 

If  we  now  expand  equation  (3.4)  a  quadratic  is  obtained  which  can  be 

written  as  ^1  ^J^"^        »  ^^^^"^   ^'[/^s^-fTl  •^^^(jtiv'-f.T))^      »   the 
initial  point  values.  Solving  this  equation  by  means  of  the  quadratic 
formula  gives 

(3.5) 

Equation  (3.5)  can  now  be  expanded  in  a  binomial  expansion  giving  in  a 
manner  similar  to  Case  It 

^^  -  S^    •  (3.6) 

The  choice  of  the  positive  sign  for  the  second  term  is  discussed  under 
Case  1. 

Equation  (3.4)  can  be  now  written  as  relative  vorticity 
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3*     Case  of  No  Horizontal  TTlnd  Shear* 

If  the  horiiontal  shear  of  the  -wind   is  either  assumed  to  be  lero  or 
is  neglected  then  the   general  equation  (2*8)   simplifies  to 

This  can  be  ^written  as 

We  see  then  that  for  the  oase  of  no  shear  the  ratio  of  the  absolute 
vortioity  at  two  points  is  equal  to  the  ratio  of  their  respective  absolute 
temperatures*     It  is  also  to  be  noted  that  this  relationship  which  involves 
no  assumptions  as  to  divergence  of  the  wind  or  barotropy  differs  from  Rossby*8 
formula  as  applied  to  the  forecasting  of  particle  trajectories  only  by  the 
presence  of  the  temperature  term* 

If  we  express  the  vorticity  as  l(^^-f-7i    aq<1  choose  an  inflection  latitude^ 
as  did  Rossby  and  co-workers     f  llj   ,   at  which  the  curvature   is   zero,  then 
equation  (3*8)  can  be  expressed  as 

^'^=^^-^  (3..) 

which  corresponds  to  Rossby* s  formula 

Equation  (3*9)  could  be  applied  in  a  manner  entirely  analogous  to  the 
technique  used  by  Rossby  in  forecasting  air  particle  trajectories.  The 
coriolls  force  at  any  point  from  the  inflection  latitude  can  be  expressed 
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as  a  function  of  distance  from  such  latitude  so  that  /)  r  A^ -f /^«  where  ^ 
is  the  rate  of  change  of  ooriolis  force  with  latitude  and  is  considered 
constant  oyer  the  latitude  range  concerned*  Then  (3«9)  can  be  written  as 

>^/'  r  2)- -  ^^4  ^  ^0  '  (3.10) 

If  we  express  the  curvature  in  terms  of  a  second  order  differential 
equation  of  its  V  coordinate  or  distance  from  the  inflection  latitude  we 
obtain  after  substituting  this  in  (3.10) 


This  differential  equation  could  perhaps  be  solved  by  elliptic  integrals 
if  a  proper  constant  is  chosen  for  the  temperature  term  so  that  a  value  of  V 
coordinate  would  be  obtained  as  a  function  of  the  intersection  angle  of 
particle  path  with  inflection  latitude • 

However,  a  qualitative  interpretation  of  equation  (3.10)  can  be  made 
showing  effects  on  the  relative  vortioity  with  varying  paths  without  actually 
solving  the  differential  equation.  If  the  particle  is  moving  northward  init- 
ially, the  effect  of  coriolis  force  is  to  decrease  the  relative  vorticity 
and  increasing  temperature  tends  also  to  increase  the  vorticity.  However, 
now  assuming  decrease  in  temperature  northward  the  temperature  term  in 
general  is  of  much  smaller  magnitude  than  f3if    ,   hence  the  particle  curves 
anticyclonically  but  to  a  slightly  greater  extent  than  the  path  resulting 
from  coriolis  considerations  only.   If  the  temperature  increases  northward 
along  the  path  the  effect  is  to  decrease  the  effect  of  increasing  ooriolis 
force  so  that  relative  vorticity  is  decreased  to  a  lesser  extent  than  when 
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the  coriolis  effect  is  considered  alone.     Therefore  the  path  would  curve 
less  anticyolonioally  than  otheirwise. 

For  a   southward  moving  current  the  converse  of  all  the  above   state- 
loents  would  hold* 
4.     Combinations  of  the  Three  Preceding  Cases* 

a*     Geostrophio  Flow  with  no  Wind  Shear* 

Under   such  conditions   of  wind  flow  equation   (2.8)  becomes  simply 

which  can  be   simplified  to 

or  the   change   of  vortioity^  which  is  now  due   simply  to  change   of 
coriolis   force  ) is  directly  proportional  to  the   change  of  respective 
absolute  temperatures*     It  is  also  interesting  to  note  that  the  flow 
in  this  case  implies  a  gradient  of  temperature  from  South  to  North, 
which  of  course,    is  not  true   in  nature  except  under   local  conditions • 
b.     Geostrophic  Flow  with  a   Barotropio  Atmosphere* 

In  this  case  the  general  equation  (2*8)   can  be  written 
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c*     A   Barotropic  Atmosphere  with  no  Vfind  Shear. 
Under  these  conditions  equation   (2*8)  becomes 


ese  oonaitions  equa 


This   can  be  written  as 
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This  result  is  identical  in  form  to  that  used  by  Rossby  and  co- 
workers [111  i^  forecasting  constant  vortioity  trajectories  with  the 
significant  difference  that  no  restrictions  have  been  placed  on  divergence 
or  convergence  in  this  development.  Therefore,  the  same  technique  as  used 
by  Rossby  in  forecasting  particle  trajectories  could  be  applied  to  the 
gradient  wind  using  identical  relationships  as: 

where  i^  s  rate  of  change  of  coriolis  force  and  (y)  equals  distance  from 
an  inflection  latitude  at  which  relative  vorticity  and  curvature  are  zero* 
Solutions  of  the  second  order  differential  equation  for  the  maximum  y  co- 
ordinate as  a  function  of  the  intersection  angle  of  the  path  with  the  in- 
flection latitude  would  be  identical  to  those  of  Rossby  f ll3  or  Fultz  C4]. 

The  verification  of  forecasts  of  such  computed  trajectories  should  be 
more  accurate  than  that  of  Rossby  because  of  the  added  refinement  of  no 
restriction  on  divergence  and  convergence* 

d.  Combination  of  all  Three  Special  Cases  —  Geostrophio  Flow,  No 
Wind  Shear,  and  Barotropy. 

These  conditions  would  give  a  form  for  the  vorticity  equation  thus 

Jl.  ?. 

which  gives  the  rather  trivial  result 

However,  absolute  vorticity  is  again  shown  to  be  constant  although  only 
a  function  of  latitude.  Furthermore,  tonal  flow  is  implied  by  this 
equality  of  coriolis  terms. 
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IV.   SOLUTION  OF  THE  VGRTICITY  EQUATION  FOR  A  CONSTANT 
PRESSURE  SURFACE 


The  equation  of  absolute  vorticity  along  a  surface  of  constant 
pressure  is  somewhat  more  simplified  than  along  a  constant  height  sur- 
face because  the  solenoid  term  is  absent*  This  is  not  surprising  since 
an  isobar ic  surface  obviously  cannot  be  intersected  by  solenoids. 

Thus,  equation  (2.2)  when  applied  to  an  isobar  ic  surface  involves 
only  a  divergence  term: 

We  can  now  express  this  divergence  in  terms  of  ooriolis  force,  cur- 
vature and  velocity  in  a  manner  entirely  analogous  to  that  for  a  constant 
height  surface* 

Starting  with  the  expression  for  the  gradient  wind: 

y=  v..  -  (^  M 

which  we  can  express  for  an  isobario  surface  as 


Now  multiplying  by  ^  and  rearranging: 
solving  for  V   gives 


(4*2) 
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Taking  the  divergence  of  both  sides  gives 
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whioh  reduces  to 

^  (4.3) 

since   Vy.\7   <J^aO  . 

It  is  seen  that  equation   (4«3)   gives  the  divergence  of  the  gradient 

wind  as  the  scalar  product  of  a  vector  parallel  to  the  isojaypses  and  the 

-/ 

gradient  of  the  quantity  (i^-f^V    . 

Upon  substituting  this  expression  for  the  divergence   in   (4«l)  and 

,  .  '\.A^.^^J     ^'•1^'y^    Uml> 

rearranging  we  get  '-  \ 

Since  •^^V^Kii-=-f-Y(^+^V)  we  can  write,  after  performing 

indicated  differential  operations: 


which  can  be  expressed  as  the  integral 


j^d(Lf^)^  7  C^^VjL(7^^i^,v)^x-, 
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Now  since  ^iH  -  V      and  Va.'  <^t-   J*<'^ 
we  can  write  (4.4)  as 

which  is  readily  integrable  to 

"*     ^   .    Q)±J^  (4.6) 
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Equation  (4.6)  is  identical  with  (3.4)  which  was  obtained  for  a 
constant  level  surface  by  assuming  a  barotropic  atmosphere*  However,  in 
developing  equation  (4.5)  no  such  assumption  was  necessary.  For  equation 
(4.5)  to  be  applicable,  however,  the  flow  must  be  assumed  to  take  place 
along  an  isobario  surface*  Therefore  its  application  would  be  highly 
limited  because  of  the  improbability  of  isobar ic  flow. 

Now  equation  (4*5)  can  be  rearranged  as  a  quadratic  and  expanded 
binomially  in  a  manner  similar  to  equation  (3*4)  and  which  gives  Identical 
results: 


f~e)  .  ^^ 


As  a  suggested  extension  to  this  work  solutions  of  the  vorticity 
equation  could  be  carried  out  in  an  entirely  analogous  manner  for  isentroplc 
surfaces  and  would  have  wider  application  than  the  solution  for  a  constant 
pressure  surface  because  of  the  greater  prevalence  of  flow  along  isentropic 
surfaces* 
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